#load packages and read in data 

library(tidyverse)
## ── Attaching packages ───────────────────────────────────────────────────────────────── tidyverse 1.2.1 ──
## ✔ ggplot2 3.2.0     ✔ purrr   0.3.3
## ✔ tibble  2.1.3     ✔ dplyr   0.8.3
## ✔ tidyr   0.8.3     ✔ stringr 1.4.0
## ✔ readr   1.3.1     ✔ forcats 0.4.0
## ── Conflicts ──────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
norep_noclim_df <- read_csv("norep_noclim.csv")
## Parsed with column specification:
## cols(
##   b = col_double(),
##   c = col_double(),
##   year = col_double(),
##   r = col_double(),
##   f = col_double(),
##   f_msy = col_double(),
##   f_ratio = col_double(),
##   f_ratio_err = col_double(),
##   id = col_double(),
##   r_0 = col_double(),
##   error = col_double(),
##   r_s = col_double()
## )
#alter dataset slightly, may be unnecessary, but we shall see 
nr_nc_df <- norep_noclim_df %>%
   mutate(parameters = paste0(paste0(paste0(paste0(paste0("r=", r), ", f="), f), ", err="), f_ratio_err)) %>% 
  group_by(id)

nr_nc_plot <- ggplot(nr_nc_df, aes(x = year, y = b, group = id)) +
  geom_line(aes(color = r)) +
  theme_bw()

nr_nc_plot

norep_clim_df <- read_csv("norep_clim.csv")
## Parsed with column specification:
## cols(
##   b = col_double(),
##   c = col_double(),
##   year = col_double(),
##   r = col_double(),
##   f = col_double(),
##   f_msy = col_double(),
##   f_ratio = col_double(),
##   f_ratio_err = col_double(),
##   id = col_double(),
##   r_0 = col_double(),
##   error = col_double(),
##   r_s = col_double()
## )
#alter dataset slightly, may be unnecessary, but we shall see 
nr_c_df <- norep_clim_df %>%
   mutate(parameters = paste0(paste0(paste0(paste0(paste0("r=", r), ", f="), f), ", err="), f_ratio_err)) %>% 
  group_by(id) %>% 
  filter(f!=0) %>% 
  filter(r_s < 0 & r_s >= -0.007) 

nr_c_df$error <- as.factor(nr_c_df$error)

nr_c_plot <- ggplot(nr_c_df, aes(x = year, y = b, group = id)) +
  geom_line(aes(color = r_0)) +
  theme_bw() +
  scale_color_continuous(low = "#079EDF", high = "#B8CE55", name = "Growth Rate") +
  labs(title = "No Adaption + Climate", y = "Biomass", x = "Year") +
  theme(axis.text.x=element_text(size=10), plot.title = element_text(hjust = 0.5, face = "bold", size = 15), axis.title.x = element_text(face = "bold"), axis.title.y = element_text(face = "bold"))

nr_c_plot

ggsave(filename="norep_clim.jpg", plot=last_plot())
## Saving 7 x 5 in image
rep5_clim_df <- read_csv("rep5_clim.csv")
## Parsed with column specification:
## cols(
##   b = col_double(),
##   c = col_double(),
##   year = col_double(),
##   r = col_double(),
##   f = col_double(),
##   f_msy = col_double(),
##   f_ratio = col_double(),
##   f_ratio_err = col_double(),
##   id = col_double(),
##   r_0 = col_double(),
##   error = col_double(),
##   r_s = col_double()
## )
#alter dataset slightly, may be unnecessary, but we shall see 
r5_c_df <- rep5_clim_df %>% 
  group_by(id) %>% 
  filter(f!=0) %>% 
  filter(r_s < 0 & r_s >= -0.007) %>% 
  filter(error == 0.3)

r5_c_df$error <- as.factor(r5_c_df$error)
r5_c_df$r_s <- as.factor(r5_c_df$r_s)

r5_c_plot <- ggplot(r5_c_df, aes(x = year, y = b, group = id)) +
  geom_line(aes(color = r_0)) +
  theme_bw() +
  scale_color_continuous(low = "#079EDF", high = "#B8CE55", name = "Growth Rate") +
  labs(title = "5 Year Adaption + Climate", y = "Biomass", x = "Year") +
  theme(axis.text.x=element_text(size=10), plot.title = element_text(hjust = 0.5, face = "bold", size = 15), axis.title.x = element_text(face = "bold"), axis.title.y = element_text(face = "bold"))

r5_c_plot

ggsave(filename="rep5_clim.jpg", plot=last_plot())
## Saving 7 x 5 in image
rep10_clim_df <- read_csv("rep10_clim.csv")
## Parsed with column specification:
## cols(
##   b = col_double(),
##   c = col_double(),
##   year = col_double(),
##   r = col_double(),
##   f = col_double(),
##   f_msy = col_double(),
##   f_ratio = col_double(),
##   f_ratio_err = col_double(),
##   id = col_double(),
##   r_0 = col_double(),
##   error = col_double(),
##   r_s = col_double()
## )
#alter dataset slightly, may be unnecessary, but we shall see 
r10_c_df <- rep10_clim_df %>% 
  group_by(id) %>% 
  filter(f!=0) %>% 
  filter(r_s < 0 & r_s >= -0.007)%>% 
  filter(error == 0.3)

r10_c_df$error <- as.factor(r10_c_df$error)

r10_c_plot <- ggplot(r10_c_df, aes(x = year, y = b, group = id)) +
  geom_line(aes(color = r_0)) +
  theme_bw() +
  scale_color_continuous(low = "#079EDF", high = "#B8CE55", name = "Growth Rate") +
  labs(title = "10 Year Adaption + Climate", y = "Biomass", x = "Year") +
  theme(axis.text.x=element_text(size=10), plot.title = element_text(hjust = 0.5, face = "bold", size = 15), axis.title.x = element_text(face = "bold"), axis.title.y = element_text(face = "bold"))


r10_c_plot

ggsave(filename="rep10_clim.jpg", plot=last_plot())
## Saving 7 x 5 in image
rep15_clim_df <- read_csv("rep15_clim.csv")
## Parsed with column specification:
## cols(
##   b = col_double(),
##   c = col_double(),
##   year = col_double(),
##   r = col_double(),
##   f = col_double(),
##   f_msy = col_double(),
##   f_ratio = col_double(),
##   f_ratio_err = col_double(),
##   id = col_double(),
##   r_0 = col_double(),
##   error = col_double(),
##   r_s = col_double()
## )
#alter dataset slightly, may be unnecessary, but we shall see 
r15_c_df <- rep15_clim_df %>% 
  group_by(id) %>% 
  filter(f!=0) %>% 
  filter(r_s < 0 & r_s >= -0.007)%>% 
  filter(error == 0.3)

r15_c_df$error <- as.factor(r15_c_df$error)

r15_c_plot <- ggplot(r15_c_df, aes(x = year, y = b, group = id)) +
  geom_line(aes(color = r_0)) +
  theme_bw() +
  scale_color_continuous(low = "#079EDF", high = "#B8CE55", name = "Growth Rate") +
  labs(title = "15 Year Adaption + Climate", y = "Biomass", x = "Year") +
  theme(axis.text.x=element_text(size=10), plot.title = element_text(hjust = 0.5, face = "bold", size = 15), axis.title.x = element_text(face = "bold"), axis.title.y = element_text(face = "bold"))


r15_c_plot

ggsave(filename="rep15_clim.jpg", plot=last_plot())
## Saving 7 x 5 in image
rep20_clim_df <- read_csv("rep20_clim.csv")
## Parsed with column specification:
## cols(
##   b = col_double(),
##   c = col_double(),
##   year = col_double(),
##   r = col_double(),
##   f = col_double(),
##   f_msy = col_double(),
##   f_ratio = col_double(),
##   f_ratio_err = col_double(),
##   id = col_double(),
##   r_0 = col_double(),
##   error = col_double(),
##   r_s = col_double()
## )
#alter dataset slightly, may be unnecessary, but we shall see 
r20_c_df <- rep20_clim_df %>% 
  group_by(id) %>% 
  filter(f!=0) %>% 
  filter(r_s < 0 & r_s >= -0.007)%>% 
  filter(error == 0.3)

r20_c_df$error <- as.factor(r20_c_df$error)

r20_c_plot <- ggplot(r20_c_df, aes(x = year, y = b, group = id)) +
  geom_line(aes(color = r_0)) +
  theme_bw() +
  scale_color_continuous(low = "#079EDF", high = "#B8CE55", name = "Growth Rate") +
  labs(title = "20 Year Adaption + Climate", y = "Biomass", x = "Year") +
  theme(axis.text.x=element_text(size=10), plot.title = element_text(hjust = 0.5, face = "bold", size = 15), axis.title.x = element_text(face = "bold"), axis.title.y = element_text(face = "bold"))


r20_c_plot

ggsave(filename="rep20_clim.jpg", plot=last_plot())
## Saving 7 x 5 in image
rep10_noclim_df <- read_csv("rep10_noclim.csv")
## Parsed with column specification:
## cols(
##   b = col_double(),
##   c = col_double(),
##   year = col_double(),
##   r = col_double(),
##   f = col_double(),
##   f_msy = col_double(),
##   f_ratio = col_double(),
##   f_ratio_err = col_double(),
##   id = col_double(),
##   r_0 = col_double(),
##   error = col_double(),
##   r_s = col_double()
## )
#alter dataset slightly, may be unnecessary, but we shall see 
r10_nc_df <- rep10_noclim_df %>%
  group_by(id) %>% 
  mutate(initial_r = ifelse(r <= .3, "slow-growing", ifelse(r >.3 & r <= .5, "medium-growing", "fast-growing"))) %>% 
  filter(f>=.1)

r10_nc_plot <- ggplot(r10_nc_df, aes(x = year, y = b, group = id)) +
  geom_line(aes(color = r)) +
  theme_bw() +
  scale_color_continuous(low = "#079EDF", high = "#B8CE55", name = "Growth Rate") +
  labs(title = "10 Year Adaption, No Climate", y = "Biomass", x = "Year") +
  theme(axis.text.x=element_text(size=10), plot.title = element_text(hjust = 0.5, face = "bold", size = 15), axis.title.x = element_text(face = "bold"), axis.title.y = element_text(face = "bold"))

r10_nc_plot

ggsave(filename="rep10_noclim.jpg", plot=last_plot())
## Saving 7 x 5 in image